Continuum generation by dark solitons 
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We demonstrate that the dark sohton trains in optical fibers with a zero 
of the group velocity dispersion can generate broad spectral distribution 
(continuum) associated with the resonant dispersive radiation emitted by 
solitons. This radiation is either enhanced or suppressed by the Raman 
scattering depending on the sign of the third order dispersion, (c) 2009 
Optical Society of America 

OCIS codes: 000.0000, 999.9999. 

One of many fundamental and practical outcomes of research into supercontinuum generation 
in optical fibers [1,2] has been a change of emphasis in studies of temporal solitons. If in 
the pre-supercontinuum era fiber solitons were mostly perceived as information carriers [3], 
now days they are widely recognized and researched as facilitators and mediators of efficient 
frequency conversion in infrared and visible parts of spectrum [1,2]. 

Trains of bright solitons naturally emerge in fibers from an intense pump pulse and their 
role in supercontinuum generation and frequency conversion has been extensively explored 
[1,2]. In particular, it has been demonstrated that the interplay of the Raman scattering 
and third order dispersion (TOD) can lead to the exponential amplification of the resonant 
(Cherenkov) radiation emitted by a soliton [4]. This happens if the central soliton frequency 
is shifted by the Raman effect towards a zero of the group velocity dispersion (GVD), which is 
possible if TOD is negative [4] . Amplification of the resonant radiation has been used to create 
coherent infrared fiber based sources [5-7]. While, the widest octave broad supercontinua 
extended from ~400nm to ~2/im have been most commonly observed across the zero GVD 
point with positive TOD [8] . Generation of these spectra is fueled by the cascaded scattering 
of the short wavelength radiation on and its trapping by the Raman solitons [2,9, 10]. 



Theory of the radiation emission by bright sohtons has been developed well before the 
recent wave of supercontinuum research [1,11]- However, importance of the Raman effect 
in the soliton-radiation interaction has been revealed and explained only recently [4,9]. 
Similarly, it is known that dark solitons also emit radiation when perturbed by TOD [12,13], 
but the interplay of this process with the Raman scattering has not been studied so far and is 
the subject of this Letter. Trains of dark solitons can be created with the established methods 
[14, 15] and, as we will demonstrate below, can generate continuum of dispersive radiation. 
Continuum generation by dark solitons can be important not only from a fundamental point 
of view, but also practically relevant if the pump frequency is in the normal GVD range of 
a fiber and generation of a broadband signal across the zero GVD wavelength is sought. 

A model we use is the dimensionless NLS equation with the TOD and Raman terms 

id,A = ^d^A + ied^A 

-(1 - e)\A\^A -OA j R{t')\A{t - t')\Mt'. (1) 

Here R is the Raman response function of silica and 9 = 0.18 [1]. t is the delayed time 
measured in the units of T = lOOfs. e = /53/[6T/52] is a small parameter characterizing the 
ratio between TOD {(3-^) and normal GVD {[32 > 0). A is the amplitude of the electric 
field scaled to Vp, where P = l/^yLn]- Lc = T'^/(32 is the dispersion length and 7 is the 
nonlinear fiber parameter. Assuming A ~ e"*''*, the zero GVD frequency is 6q = — l/[6e]. 
For positive (or negative) TOD the anomalous GVD range hostile to the existence of dark 
solitons is located for 6 < 60 {or 6 > 60). A unit of 6 corresponds to lOTHz. 

We are interested in A consisting of the dark soliton F plus small radiation g: A{z,t) = 
{F{t, z) + (^(r, z))e^'^^^^ where t = t — zKsincj). The soliton solution is [16] 

F = K tanh[rKcos0] COS0 — msin0. (2) 

K > is the amplitude of the soliton background. is the darkness parameter, with the 
condition cos0 > implying G [— 7r/2, 7r/2]. = corresponds to the black soliton 
having zero amplitude in the middle and propagating with the group velocity matching the 
background velocity [16]. G (0,7r/2) and G (— 7r/2,0) correspond to the dark solitons 
propagating, respectively, slower and faster than the background wave. 
Solutions of a linearized equation for g are sought in the form 

Note, that for A ~ e"*''*, the exponent e*'^'''^ with 6r > corresponds to the overall negative 
frequency shift. Condition A = leads to the resonant frequency [12, 13] 

6j. = — ^ [a + v^aM-^K^e^75os^] , a = 2eK sin + -. (4) 
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Expanding in e we find 5r = l/(2e)(l + 4Kesin0 + 2K^e^cos^0 + 0(e^)), so that l^^l > l^ol, 
i.e. the resonant frequency is detuned further away from the sohton, than the zero GVD 
frequency. An important difference with bright solitons is that the amplitudes Gi and G2 
are not independent here, but are coupled by nonlinear interaction mediated by the finite 
amplitude soliton background. The ratio of these amplitudes works out as 

Gi ^ 1 

Eq. (3) implies that the radiation spectrum is expected to peak on both sides from the 
soliton spectral center at 5 = 0. Since \5r\ > \So\, one of the radiation peaks belongs to 
the normal GVD range and the other one to the anomalous GVD. From Eq. (5) it follows 
that the amplitude of the radiation belonging to the anomalous GVD range is much larger. 
The above holds for both signs of TOD, i.e. e. Thus the radiation is primarily emitted 
into a spectral range where the balance of GVD and nonlinearity is not able to support 
dark solitons. Figs. l(a,b) show solitons and their radiation using XFROG spectrograms [1] 
calculated with the Raman effect disregarded (6* = 0). For e > (Fig. 1(a)) and e < 
(Fig. 1(b)) the strong and weak radiation peaks are swapped, so that the strong one always 
belongs to the anomalous GVD range. Spectrograms used here are computed as S{6, t) = 
I J A{t')sech.{t' — t) exp[—i6t']dt'\'^ . For kcos</) = 1 and t = the dark soliton spectrogram, 
A = F, can be calculated explicitly 

S (X (5^ cos^ + 25 sin cos + sin^ 0) sech^[7r5/2]. (6) 

Taking account of the Raman scattering significantly changes the radiation pattern, cf. 
Fig. 1 and Fig. 2. The Raman effect enhances both radiation peaks for positive TOD and 
suppresses them for negative TOD. Physics behind this observation is as follows. Spectrogram 
S{6, 0) of the ideal black, = 0, soliton consists of the two symmetric in 6 peaks (not 
clearly seen in Fig. 1, but obvious from Eq. (6)). It means that the resonant radiation in 
the anomalous GVD range is excited by the black soliton with equal efficiency for either 
positive or negative TOD. The Raman scattering however transfers more energy into the 
lower frequency [5 < 0) side of the soliton spectrum. It makes more of the soliton energy to 
flow into the anomalous GVD range, providing TOD is positive. For negative TOD, however, 
the effective spectral center of mass of the soliton is shifted by the Raman effect further into 
the normal GVD range, where little disturbs the balance of normal GVD and nonlinearity. 
Naturally, this depletes the energy transfer into radiation. Asymmetry of the dark soliton 
spectra induced by the Raman effect assumed in the above discussion can be readily verified 
numerically and was experimentally measured 20 years ago [17]. 

The soliton spectral asymmetry induced by the Raman scattering can be mimicked by 
taking the fast dark solitons (0 < 0), see Eq. (6). Equally it is true to say, that the Raman 
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scattering accelerates dark solitons by inducing the adiabatic drift of in the direction of 
— 7r/2. Thus by taking initially a faster (slower) soliton in a positive (negative) TOD fiber 
one shifts more of the soliton energy to the anomalous GVD range, thereby boosting energy 
transfer into dispersive radiation. This point is further discussed and illustrated below. 

A train of dark solitons is an obvious candidate for generation of multiple radiation peaks 
forming a continuous spectrum. One of the straightforward ways to generate these trains 
is through the dispersion induced interference of two suitably delayed short pulses [14,15]. 
Each interference dip has its own carrier frequency associated with the chirp developing 
across the whole interference pattern. The characteristic GVD and nonlinear lengths can be 
made approximately equal and are typically shorter then the TOD length. Thus one can 
arrange initial conditions so that the solitons are formed first and then the TOD induced 
radiation starts to develop. 

Fig. 3 shows time domain evolution of the soliton train for e > (a) and for e < 
(b). Both are with the Raman effect included. The solitons are separated into two groups. 
The ones found for t < and t > are the fast (0 < 0) and slow (</> > 0) solitons, 
respectively. For e > the spectra of the fast solitons are shifted towards the range of 
anomalous GVD. Therefore these solitons are more keen to emit radiation, than the slow 
solitons (see pronounced oscillations developing for t < in Fig. 3(a)). The Raman effect 
enhances the soliton spectral asymmetry further and boosts the radiation emission by the 
fast solitons, see Fig. 3(b). For e < spectra of the slow solitons have greater overlap with 
the anomalous GVD range. However, the Raman effects tends to compensate for this. Thus 
the net effect is the radiation depletion, cf. Fig. 3(a) and (b). 

Spectral evolution along the fiber length for the cases of positive and negative TOD cor- 
responding to the time plots in Fig. 3 are shown in Fig. 4. Far stronger energy transfer from 
the normal to the anomalous GVD range is obvious in the positive TOD case. Quantitative 
characterization of this effect is facilitated by Fig. 5, where the full lines show the spectral 
distributions (plotted using the log-scale) corresponding to the final time domain signals in 
Fig. 3. The spectrum generated in the anomalous GVD range is few orders of magnitude 
stronger, than the one in the negative TOD case. Contrary, if the Raman effect is switched 
off, then e > case is a mirror image of the e < case, cf. the dashed lines in Fig. 5(a) and 
(b). 

In summary, we have demonstrated that the Raman scattering significantly amplifies res- 
onant dispersive radiation emitted by a train of dark solitons in fibers with the positive third 
order dispersion. This leads to generation of a broad band continuum in the anomalous GVD 
range, when a pair of pump pulses is applied in the normal GVD. Practical realizations of 
this effect and its applications for frequency conversion should be investigated further. 
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Figure captions 

Fig. 1. (Color online) XFROG spectrograms showing resonant radiation by a dark soliton 
without the Raman effect = 0). Dashed vertical lines indicate ±5^ predicted by Eq. (4). 
Full vertical lines indicate the zero GVD frequency 6q. (a) e = 0.0833 and (b) e = —0.0833. 
Other parameters are k = 1, (p = 0, z = 30. 



Fig. 2. (Color online) The same as Fig. 1, but with the Raman effect {6 = 0.18). Strong 
emission of radiation by the black soliton, see (a), is accompanied by creation of a shallow 
dark soliton with (/) close — 7r/2 (marked with an arrow). 



Fig. 3. (a) Time domain evolution of the dark soliton train in a fiber with positive TOD: 
e = 0.0217, 9 = 0.18. Initial condition is A = VTO[sech{t - 3) + sech{t + 3)]. (b) is the same 
as (a), but for negative TOD: e = -0.0217. 



Fig. 4. (Color Online) Spectral development of continuum generation by a train of dark 
solitons corresponding to the time domain evolution in Fig. 3, but for larger z. Full vertical 
lines indicate the zero GVD frequency. 

Fig. 5. Full lines in (a) and (b) show spectra corresponding to the 2; = 10 time domain signals 
in Fig. 3. Dashed lines show the corresponding spectra computed without the Raman effect 
(9 = 0). 
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